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Massless Dirac particles cannot be confined by an electrostatic potential. This is a problem for 
making graphene quantum dots but confinement can be achieved with a magnetic field and here, 
general conditions for confined and deconfined states are derived. There is a class of potentials 
for which the character of the state can be controlled at will. Then a confinement-deconfinement 
transition occurs which allows the Klein paradox to be probed experimentally in graphene dots. 
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Single layer graphene is attracting attention because 
its charge carriers are massless, relativistic particles [ij. 
The relativistic effects result from a unique, zero-gap 
band structure that leads to quantum states described 
by the two-component Dirac- Weyl equation. This allows 
relativistic physics to be explored in a solid state system 
and has many potential applications ranging from high 
frequency electronics f]\ to quantum computing [2]. In 
particular, graphene quantum dots are very attractive as 
spin qubits because they are expected to have low spin 
decoherence 0. However there is a problem with mak- 
ing graphene dots for quantum computing or any other 
application: relativistic effects prevent massless particles 
from being confined by an external potential. 

The problem results from the Klein paradox [1, U, @ . 
When relativistic particles with mass are incident on a ID 
potential barrier, the state in the barrier decays exponen- 
tially unless the barrier height exceeds the threshold for 
pair production, at which point the state in the barrier 
becomes oscillatory. The paradox is that any attempt to 
enhance the localisation by increasing the barrier height 
eventually destroys it. But there is no threshold for pair 
production for massless particles so exponential decay 
and bound states do not occur. 

Graphene dots can be formed from external poten- 
tials or nanocrystals. The quantum states [1, 0, B 
in external potentials are quasi-bound: they have 
a low amplitude oscillatory tail and are similar to the 
scattering resonances studied in undergraduate physics. 
A perpendicular magnetic field enhances the localisation 
of these states fi\ and true bound states can occur in 
graphene dots defined by a spatially non-uniform field 
3 . So a magnetic vector potential has a localising effect 
that tends to cancel the delocalising effect of a scalar po- 
tential. But what are the general conditions for confined 
states to occur when an electrostatic scalar potential and 
a magnetic vector potential are applied to graphene si- 
multaneously? 

This question is answered in the present Letter. It 
is shown that both true bound states and quasi-bound 
states occur, depending on the form of the potentials. In 
addition, there is a third and most interesting possibil- 
ity. In some cases, the character of the states depends on 



the parameters of the potentials and can be controlled 
at will. A confinement-deconfinement transition then oc- 
curs in which the character of the states changes from 
oscillatory to exponential as in the Klein paradox for par- 
ticles with mass. This gives a way of probing the Klein 
paradox experimentally in a solid state system and nu- 
merical studies of the quantum states in a realistic dot 
model show it is feasible. 

The system considered here is a cylindrically symmet- 
ric graphene dot. The dot region is defined by an electro- 
static potential, V{r), and the magnetic vector potential, 
Ag(r), is in the azimuthal (9) direction. The functional 
form of these potentials has a significant effect on the 
character of the states. In particular, when V and Ag 
increase as power laws, V = Vbr*, Ag ~ Aqt^, s,t > 0, 
the character of the states depends critically on s and 
t. li s > t the states oscillate in the asymptotic regime 
of large r but decay exponentially when s < t. In both 
cases the asymptotic character of the states is indepen- 
dent of Vq and Aq but when s = t the character of the 
states does depend on these constants and a transition 
from exponential to oscillatory behaviour occurs when Vq 
is increased or is decreased. This is the confinement- 
deconfinement transition. 

The two-component envelope function, tp, satisfies 
{V + jcr ■ 7r/h)ip — Eip, where cr are the Pauli ma- 
trices, TT — p + eA and 7 = 646 meV nm [7]. For 
cylindrically symmetric systems ^ = (xi (?') exp(i(m — 
1)0), x^if) exp{im9)) where m is the angular momentum 
quantum number. The radial functions xi ^-nd X2 satisfy 
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These equations can be uncoupled by differentiating 
them and this leads to the relation 



X2 + a{r)x2 + bir)X2 = 0, 



where 



1 1 dV 



(3) 



(4) 



h dr 7^ fi'^ 

The first order derivative in Eq. ([3]) is eliminated by 
putting X2ir) — U2{r) exp(— j a{r)dr/2) which gives 

^2+^2(0^^2=0, (6) 

where fc|(r) — b—a'/2 — a?lA. Although fcf diverges when 
E X2 is regular there. Since exp(— / a{r)dr/2) is 

not an oscillatory function of r, has the same charac- 
ter as X2 • Eq. ^ shows that this character is oscillatory 
when fc| is positive and exponential when fc| is nega- 
tive. Asymptotic exponential decay is characteristic of 
a bound state but here "confined" is used to indicate a 
bound state that is localised near the centre of the dot. 

When V and Aq increase as power laws, the asymp- 
totic form of fc| is (Vo/7)^r^'' — (e^o/^)^'^^*- Hence A:| is 
positive when s > t, leading to oscillatory character and 
deconfined states. And /c| is negative when s < t, leading 
to exponential character and confined states. However if 
s = t the asymptotic form is [(Vo/7)^ — {eAo / h)'^]r'^* so 
the sign of A:| depends on Vq and Aq and a confinement- 
deconfinement transition occurs when Vq = (76^0/^)^- 
This also follows from the uncoupled equation for xi- 
the corresponding /c^-value, ^ fc| but kf A;| in the 
asymptotic limit so xi E^nd X2 have the same character. 

To investigate the states further, Eqs. ^ and ^ are 
solved numerically. The Hamiltonian, H , satisfies 

-271-^7 [(xi„X2/3 + xLxi/3)r]^ , (7) 

where ijja and ipp are two-component states of angular 
momentum m and Xia^Xi|3 ^-re the corresponding radial 
functions. Eqs. ([1]) and ^ lead to a Hermitian eigenvalue 
problem when the boundary terms in Eq. ([7]) vanish. For 
this to happen it is sufficient that one component is reg- 
ular at the origin and one component vanishes at the 
boundary, R. Then it follows from Eqs. H]) and ^ that 
both components are regular at the origin but it does not 
follow that both components vanish at R. However a true 
bound state has an exponential tail so both components 
of a bound state vanish in the limit of large R. Hence 
confined and deconfined states can be distinguished by 
solving Eqs. ([T]) and ^ subject to the boundary con- 
ditions that one component is regular at the origin and 
one component vanishes at R and then looking for an 
exponential tail in both components. 

Eqs. ^ and ^ are solved by discretizing them on 
a uniform grid. By applying the time reversal operator 
to Eqs. (H]) and ([2]) it can be shown that £^(m,Ao) = 
E{\ — m, — Aq). It is important to ensure that the eigen- 
values of the discretized Hamiltonian have the same prop- 
erty and this requires identical numbers of grid points for 
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FIG. 1: Confined and deconfined states for cases when no 
transition occurs. The frame above each state shows fcf . Top- 
most frame: E{B) for s — 2, t = 1. 



Yi and X2- This excludes the use of centred differences 
[ll| so d/dr is approximated by the forward difference op- 
erator Lf or the backward difference operator Lf,. The 
procedure depends on m. When to < 0, X2(-R) is cho- 
sen to vanish [l^, Lf is used to find dx2/dr and Lb to 
find dxi/dr and vice- versa when m > 1. Although this 
guarantees that Eim^Ao) = E{1 — to, — Aq), it has the 
disadvantage that numerical errors are linear in the step 
length Ar. To compensate for this, Ar is kept small 
and all the eigenvalues computed in this work are accu- 
rate to about 1 part in 10'^, except for s = 2,i = 1 and 
s = 2, t = 2, where there are rapid oscillations but the ac- 
curacy is still better than 2%. The discretization leads to 
a non- Hermitian matrix eigenvalue problem. A similarity 
transformation is used to reduce this to a real, symmetric 
eigenvalue problem which is solved numerically. 

The quantum states shown in the present work are all 
rn — 1 states. Other states exhibit similar features, al- 
though the amplitude of oscillations is m-dependent 
Since the main focus of this work is on the confinement- 
deconfinement transition, the states are selected so that 
the behaviour at the origin does not change significantly 
when the potential parameters are changed. All the 
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FIG. 2; As Fig. [T] but for two cases which exhibit a 
confinement-deconfinement transition. 



states have been selected to have a large amplitude close 
to the origin and can be regarded as dot states. The ener- 
gies as a function of magnetic field typically show a series 
of anti-crossings (Fig. [T|) and where necessary, the line of 
anti-crossings is followed to preserve the qualitative form 
of the state at the origin ^13i | . Every energy computed 
here is between 26 and 67 meV, within the validity limit 
of the linear graphcne Hamiltonian (« ±1 eV 

Confined and deconfincd states are illustrated in Fig.[T] 
The radial probability distribution, = r|xip, i — 1,2 
is shown together with kf. The insets show on a 
logarithmic scale. R — 600 nm, large enough to en- 
sure that the asymptotic sign of fc| has been reached. 
The magnetic field, B, is uniform. When s = 2, t = 1, 
Vq — 5 X 10^^ meV nm^^, the asymptotic sign of is 
positive and the asymptotic character of j/^p is oscilla- 
tory independent of B. The amplitude of the oscillations 
decreases with increasing B. When s — 0.5, t = I, Vq = 5 
meV nm"-'^/^ and B ^ 0, the asymptotic sign of kf is neg- 
ative and the asymptotic character of is exponential, 
as can be seen from the insets. 

The confinement-deconfinement transition is illus- 
trated in Fig.m R — 600 nm, again large enough to reach 
the asymptotic regime. When s = l,t=l,Vo = 0.5 meV 
nm~^ and B = 0.3 T, the asymptotic sign of kf is positive 
and the asymptotic character is oscillatory. In contrast, 
when B = 1.2 T, the asymptotic sign of kf is negative 
and the asymptotic character is exponential. The tran- 
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FIG. 3: V{r) for the gate geometry shown in the bottom inset. 



sition also occurs in non-uniform magnetic fields and it 
may be possible to generate a suitable field by putting a 
dot under a superconducting obstacle in a uniform field. 
Fig. [2] (bottom) shows the transition for s = 2, t = 2, 
Vq = 5 X 10~^ mcV nm^^, that is parabolic confinement 
in a linearly increasing field, B{r) =^ B^r. 

In any real quantum dot, the scalar potential would 
approach a finite asymptotic value instead of increasing 
without limit. Consequently, all the states of a graphene 
dot in a magnetic field have an exponential tail. How- 
ever, an effect similar to the confinement-deconfinement 
transition occurs in the middle distance region between 
the centre of the dot and the asymptotic regime. 

For this transition to be observable, the dot level has to 
be in the region of very low density of states between the 
bulk Landau levels. This requires a potential similar to 
the one shown in Fig. [S] The asymptotic value of the po- 
tential is engineered to be just below the dot level. This 
puts the dot level between the zeroth and first Landau 
levels. Thus the dot level can be isolated from the bulk 
Landau levels provided that they are narrow enough. 

The required potential can be generated by gate elec- 
trodes. One possible arrangement is a metal plate with 
a circular hole that contains an electrode. The graphene 
sheet is above these electrodes on 300 nm of Si02 on a 
back-gate at V. The plate (-1 V) generates the asymp- 
totically fiat part of the potential, the hole generates the 
barrier and the central electrode (-2 V) generates the 
well. Similar gated, monolayer graphene nanostructures 
have been fabricated recently pJi]. The potential in Fig. [3] 
was computed by solving the Poisson equation on a dis- 
crete mesh. Screening by the graphene sheet was treated 
in the Thomas- Fermi approximation |15i] . The resulting 
potential is magnetic field dependent because the den- 
sity of states is field dependent. This causes steps in the 
potential (inset to Fig. [3]) which occur when the num- 
ber of occupied Landau levels changes. However the field 
dependence is weak at the low fields considered here. 

The "confinement-deconfinement" transition for the 
potential in Fig. [3] and a uniform magnetic field is shown 
in Fig. [4] For all fields, the states have a peak near the 
centre of the dot and an exponential tail. The transition 
occurs in the middle distance region between these two 
features, 200 ^ ?' ^ 400 nm. For any potential, this re- 
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FIG. 4: Confinement-deconfinement transition. AsFig.[2]but 
for the potential shown in Fig. O 



gion can be identified by computing kf. In Fig. [4l this 
varies rapidly because of the steps in the potential but 
remains positive in the middle distance region when B ~ 
0.05 T. At this field, the oscillations in correspond to 
those in Fig. [2] but are much less rapid because V, hence 
kf is smaller. As the field increases the region of positive 
kf shrinks and a transition to exponential behaviour oc- 
curs. This is analogous to the change of character seen 
in the Klein paradox for relativistic particles with mass. 

The occurrence of the transition is insensitive to the 
electrode geometry. The one in Fig. [3] has the advantage 



that the graphene is easy to access but may be difficult to 
fabricate. However, similar transitions occur in systems 
with disk or spherical central electrodes with the plate 
and central electrode either above or below the graphene 
sheet. In addition, calculations for model potentials with 
a well, barrier and flat portion show that the occurrence 
of the transition is insensitive to the model parameters. 
The only requirement is a region where kf changes sign 
when ^0 increases. This is relatively easy to arrange. 

l/ip in Fig. [4] decreases by 3-4 orders of magnitude at 
r « 300 nm when the character of the state changes from 
oscillatory to exponential. This large effect could be used 
to probe the transition experimentally. For example, the 
decrease in causes a decrease in the local density 
of states (LDOS) near the dot which could be detected 
with scanning tunnelling microscopy. The decrease in 
I Zip would also cause a large decrease in the overlap of 
the dot state with states in contacts at r « 300 nm. This 
could be detected by looking for a change in transport 
through the dot state via diametrically opposite contacts 
at r « 300 nm. Numerical calculations of the LDOS for 
the potential in Fig. [3] show that the state in Fig. 0] lies 
in a region of very low bulk density of states, with the 
dot state around 0.1-0.2 meV away from any other levels. 
This suggests the dot state can be resolved experimen- 
tally in graphene samples of sufficient quality. 

In summary, confined and deconfined states occur in 
graphene dots in a magnetic field. The states are con- 
fined when the scalar potential rises slowly and decon- 
fined when it rises rapidly. But when the scalar and 
vector potentials have the same power law, the character 
of the states can be controlled at will and a confinement- 
deconfinement transition occurs. A similar effect occurs 
in a realistic dot model. This is analogous to the relativis- 
tic Klein paradox and may be experimentally observable. 
The proposed system may be attractive because it allows 
graphene dots to be formed in a uniform magnetic field. 
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